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Abstract 

We study chaotic properties of eigenstates depending on the degree of complexity in bound- 
aries of a 2D periodic billiard. Main attention is paid to the situation when the motion of a 
classical particle is strongly chaotic. Our approach allows to explore the transition from deter- 
ministic to disordered chaos, and to link chaos to the degree of ergodicity in eigenstates of the 
billiard. We have found that bouncing balls strongly reduce chaotic properties of eigenstates, 
thus leading to a serious problem in statistical description for global properties of eigenstates. 
A quite unexpected effect of rough surfaces on the form of eigenstates has been discovered and 
Oh! explained by a strong localization of a subset of eigenstates in the energy representation. 
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> ■ 1 Introduction 

in 

Electromagnetic, acoustic, or quantum waves are essential means for communications and mea- 
' surements and hence their propagation and scattering properties are of great scientific as well as 

. technological interest. During the last decade or so, there has been a surge of theoretical and exper- 

' imental breakthroughs inspired mostly by the phenomenon of localization in disordered media [1, 2]. 

. Localization is a wave interference effect and hence there is a host of potentially rich phenomena 

including, weak localization of light, Anderson localization and transport of electromagnetic radia- 
^ . tion in random media, localization of surface gravity waves on a random bottom, and electron/wave 

Vh ' propagation in films or waveguides with rough surfaces [3, 4, 5]. The well known result from the 

. scaling theory of localization [6] is that for two dimensional samples of infinite extent, there will be 

' no propagation of waves and the stationary solutions are exponentially localized if the disorder is 

O ■ strong enough. For one dimension, any amount of disorder is sufRcient. More recently there have 

' been studies on correlated random models. In one dimension, the latter allows for the appearance of 

the mobility edge, envisioned more than a decade ago but just theoretically [7] and experimentally 
[8] established a few years ago. Disorder can be caused by the medium or by the boundaries between 
media as may be the case for ballistic mesoscopic systems or for gravity waves at the surface of the 
' water in the oceans. 

Localization is known to occur not only in random models but also in deterministically quantum 
chaotic systems [9, 10]. This type of localization is known as dynamical localization since it is entirely 
due to the deterministic chaotic dynamics, where random or disordered potentials are absent. In the 
regime of fully chaotic dynamics (hard chaos) some properties of these systems may be described 
rather well by the same mathematical methods of disordered systems. In particular, by exploiting 
the nonlinear sigma model, Efetov [13] has shown that on long time scales such properties as level 
statistics of weakly disordered systems are universal and coincide with those of Random Matrix 
Theory (RMT). The latter are known to serve for description of quantum systems that are classically 
chaotic [10, 11, 12]. Additionally, based on Gutzwiller's formula [14], semiclassical calculations 
[15] predict new specific correlations in energy spectra, similar to what was obtained by another 
mathematical tool from disordered systems, the diagrammatic perturbation theory [16]. 

While there are some properties that are shared by disordered and chaotic systems, the common 
approaches employed are very good insofar as they reveal universal properties under certain condi- 
tions. Thus, theories such as RMT or the nonlinear sigma model are unable to describe particular 
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Figure 1: Geometry of the modulated billiard for Nt = 100. The inset shows one period of the 
bilUard for Nt = 6. 

features (or deviations from the universal predictions) that will distinguish between one chaotic sys- 
tem and another. One important advantage of a deterministic approach to quantum chaos is that 
it gives us the possibility of understanding particular quantum features in terms of its underlying 
classical dynamics. It is well known, for example, that scars and bouncing ball states [17] are re- 
sponsible for deviations from the predictions of Schnirelman's theory [18] for wave functions and 
also from the level statistics of RMT [19]. Further, in [20, 21, 22] it is shown that specific features 
of the structure of wave functions of quantum chaotic system are manifestations of certain type of 
periodic orbits of the underlying classical dynamics. 

In this work we examine the transition from chaos to disorder in a completely deterministic 
manner. We study the solutions of the Scliroclinger and Newton equations of a system which is 
deterministically chaotic in the classical limit, then we examine the transition chaos-disorder by 
gradually increasing the degree of complexity of the "interaction" potential till the scattering po- 
tential can be considered practically as a random one. As a paradigm we use a 2D chaotic billiard 
that is periodic in one direction. Specifically, the billiard is defined by two parallel hard walls, a flat 
one at y = and a modulated one given hy y = d + a£,{x) , where ^ is periodic ^{x + b) = ^(x) in the 
longitudinal coordinate x with period b, d is the average width of the billiard and a is the modulation 
amplitude. Here x, y, d, a, and b are dimensionless quantities obtained from Y = D + A^{2'itX/ B) 
scaled to the period B. Then, x = 2'kX/B, y = 2t:Y/B, d = 2ttD/B, a = 2'kA/B, and b = 2-k, 
where X , Y, D, A, and B have dimensions of length. The case ^(x) = cos(x), the so-called "cosine" 
or "ripple" bilhard [20, 21, 22, 23, 24, 25, 26, 27, 28], is known to yield the generic transition to chaos 
of Hamiltonian systems. Various transport and dynamical properties of the cosine billiard have been 
investigated in both, classical [23, 27, 28] and quantum regimes [20, 21, 22, 24, 25, 26, 27, 28] in its 
two versions: finite [23, 25, 26, 27, 28] (one or few periods of modulation) and infinite [20, 21, 22]. 
Here we are interested in studying the properties of the eigenstates as a function of the complexity 
of the top profile. For this purpose we write 



where A^ are random numbers in the range [-1,1]. This form allows us to choose the desire degree 
of complexity of the top boundary; from smooth {Nt 1) to extremely rough {Nt ~ 100), see 
Fig. 1. We remark that even for a single harmonic [iVr = 1 in (1)] the dynamics can yield mixed 
or full chaos, depending on the values of the parameters a, and d. Here we consider only the case 
of full chaos. Note also that | d^{x)/dx \ increases in magnitude as Nt increases, unlike the rough 
profile used in [29] where the amplitude of each harmonic N in the sum (1) is inverse proportional 
to N. Therefore, in our case for any given Nt the fluctuations are sharper, more dense, and more 
"random" than that used in [29] . 

This paper is organized as follows. In the next section we briefly review the formalism that enables 
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us to incorporate into the Hamiltonian, as an effective potential, all scattering effects that are due 
to the modulated boundary. There, we also examine the structure of the Hamiltonian matrix, from 
which we can anticipate the existence of extremely energy-localized states. In Section 3, we present 
our calculations of the energy-localization lengths for eigenstates and study them as the complexity 
of the profile increases. In Section 4, wc explore the structure of the energy-localized states in the 
case of Nt = 100, when the top profile may be considered as random, or rough. In contrast, in 
Section 5 we discuss the structure of the energy- extended states as the profiles changes from smooth 
to rough. In Section 6 we study how the energy-localized and extended states manifest themselves in 
the energy-band structure. Finally, in Section 7 we give some concluding remarks. 



2 Structure of the Hamiltonian Matrix and Bouncing Ball 
States 

Much information can be extracted from the knowledge of the structure of the Hamiltonian matrix 
which, of course, depends on the choice of basis used to expand the states. Thus, it is necessary to 
explain briefly the procedure used in the construction of the Hamiltonian matrix. 



2.1 Coordinate Transformation and Hamiltonian Matrix 

We first perform a transformation to new coordinates, {x, y) — > {u, v), such that all effects due to the 
modulated boundary become incorporated into the Hamiltonian operator. This can be accomplished 
by the curvilinear coordinates u = x and v = yd/[d + a^{x)] = y/[l + e^{x)] (e = a/d). In these 
coordinates the boundary conditions become ip{u,v) = at = 0, 1 and the Hamiltonian attains 
the form H = {l/2me)g~^/^Pag"^g^^'^Pi3g~^^^ {c(,P = u,v), where g""^ is the metric tensor given 

by 



1 + 

1 + {ve^u? 



(2) 



Here ^„ = d^/du and g = (1 + e^)^ is the metric, see Ref. [24] for details. Further, H can be 
splitted into an unperturbed part H° = (l/2me)(-P^ + Py) and an effective interaction potential 
eV{u,v,Pu,Pv)- 

Since the Hamiltonian is periodic in u, the energy eigenstates are Bloch states: 'iI)e{u,v) = 
cx.p{iku)ibkiu, v), ipk{u, v) = ijjk{u + 6, f ), where k is the quasi-momentum chosen to lie in the first 
Brillouin zone —1/2 < k < 1/2. The a-th eigenstate of energy E°'{k) can be expanded as 
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where 



^{u,v) = < u,v \ m,n >k= 
are the eigenstates of with eigenvalues 
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The resulting Hamiltonian matrix in the basis of is [24] 
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Figure 2: Central part of the 4030 x 4030 Hamiltonian matrix Hi^ii{k) with Nmax = 32, and 
Mjnax = 62. The 62 x 62 blocks corresponding to n, n' = [0, 10] are shown for the cases (a) Nt = 1, 
(b) Nt = 2, (c) Nt = 4, and (d) Nt = 6. 
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In a numerical study the Hamiltonian matrix is truncated: — A^max ^ n < N^^x^ I < m < M^^x- 
The structure or form of -ff^„ ^/„/ =< m', n'\H\m, n > depends on the ordering of the unperturbed 
basis \l >= |TO,n >. That is, wc must assign an index labeling the basis state |/ >fe= \m^n >k, 
to each pair (m, n) of indices. The spectra is, of course, independent of the assignment but there 
are two different assignments which are more convenient for our analysis. To study the structure 
of the eigcnstatcs, the most physical assignment (m, n) I is, as it will become clear below, such 

that the unperturbed energy spectrum is "energy ordered": i^j^i > -^fi ' = ^i^, On the 

other hand, to analyze the properties of the Hamiltonian matrix the more convenient assignment 
is as follows. Starting from the lowest value of n [n = ^Nmax) wc run through all values of m 
(1 < TO < Mmax)- This givcs I = 1,2, ... ,Mmax- Then we do the same for n = —N^ax + 1, 
which gives M^ax + 1 < ^ < ^Mmax, and so on untill we reach n = N^ax to get 1 < Z < L^ax, 
where L^ax — i'^Nmax + ^)Mmax is the size of the Hamiltonian matrix. This rule results in a block 
structure of Hi^i> containing {2Nmax + 1) x {2Nmax + 1) blocks, labeled by the pair (n, n') and each 

of size Mmax X Mmax- 

Fig. 2 shows the central part of the 4030 x 4030 Hamiltonian matrix [{Nmax^ M„iax) = (32,62)] 
for the cases Nt = 1, 2, 4, and 6 harmonic terms in the profile. A quick look at these figures reveals 
that the size of the bandwidth increases as Nt increases. By considering the matrix elements Hw 
of zero and first order in e [see Eq. (6)], it can be shown that the band size width AW increases 
with Nt as 

AW = {2Nt + l)Mmax, e < 1. (8) 

Moreover, the matrix elements with coefficients give rise to a secondary band that is twice as 
wide as AW, but whose elements lie mainly along the diagonal of the blocks defined by \n' — n\ = 
2N, N = [1, A'^r]. Inspection of Figs. 2(a-d) shows that such estimates agree very well with the 
numerically computed Hamiltonian matrices. The analysis also shows that the secondary band 
becomes more important as the number of harmonics Nt increases. This effect is mainly due to the 
increase, as Nt increases, of the derivative term appearing in the integrand of J^, „ [Eq. (7)] 
which comes with all terms in Hi^f. If the amplitudes An, defining the profile (1) were inversely 
proportional to TV (as in Rcf. [29]), then the corresponding Hamiltonian would couple roughly half 
as many unperturbed states. Thus, the degree of complexity of the boundary profile is manifested 
in the number of unperturbed states required to construct the eigenstates of the full Hamiltonian 
H = H° + eV. 
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2.2 Bouncing Ball States 

Another very important dynamical feature, generic of this class of billiards with small amplitude 
deformations, can be deduced from the structure of -ff^n m'n' ^"7 observing that the matrix elements 
< n,m = l\H\n' , m' = 1 > S„n' couple very weakly to the rest of the unperturbed states [sec white 
lines delimiting the blocks in Fig. 2(a-d)]. This fact results in a set of eigenstates oi H = H° + eV 
that are practically the same as the eigenstates of [see Eq. (4)] (unperturbed states) with m = 1. 
We can understand the physical origin of these weakly perturbed states, which we shall refer to as 
m = 1 states, as follows. An unperturbed eigenstate (that of a flat channel) has a "transversal 
wavelength" A^^ = 2-K/ky = 2-^ / {mn / D) = 2D/m. Then for m = 1 states, Xy is so much larger 
than the amplitude of the ripple (Xy/a. = 2/e 3> 1) that the state docs not resolve the ripplcncss 
and hence it is hardly perturbed. On the contrary, unperturbed states with large values of m are 
strongly affected, and many unperturbed states participate in the formation of such eigenstates of 
H. These expectations were numerically tested for the case Nt = 1 in Refs. [20, 21, 22]. It was 
shown that in any given range of energy one can find m = 1 states amidst others which are extended 
throughout the whole energy shell. 

Since the m = 1 states exist at all energies, it is interesting to determine the fraction of these 
states relative to the total number of states N{E) as a function of energy. This question is important 
since in view of the Schnirelman theorem [18] stating that most eigenfunctions of a classically chaotic 
(ergodic) system are equidistributed over the energy shell as ?i — > 0. In other words, we have 

lim ^-rif =0, (9) 

E^oo N{E) ' ^ ' 

where Nm=i{E) is the number of m = 1 states up to some energy E. Neither the theorem nor Eq. 
(9) says anything about how the limit is approached. Moreover, the Schnirelman theorem does not 
consider the possible existence of parabolic fixed points. These occur in our billiard; they are the 
bouncing ball orbits, or the continuous set of all horizontal trajectories {Py = 0) which do not hit 
the boundaries. Quantum-mechanically, there are no Pj^ = states but the eigenstates of the flat 
channel with minimum value of Py are the m = 1 states. Thus we sometimes refer to the m = 1 
states as "bouncing ball states". For e <C 1 one can use the unperturbed spectra given by Eq. 
(5), to get a good estimate for Nm=i{E) and N{E). For large Nmax and Mmax the relation (5) 
represents, ignoring the fc^ term, half of an ellipse in the n — m plane (recall that m G [1, M^nax] and 
n E [—N„iax, Njnax])- Thcn for large E, the number N{E) equals the area of the half ellipse, and 
the number Nm=i{E) of bouncing ball states equals twice the size or its minor axis: 



IF) 



Nm=l={-^] (11) 



therefore, 



Nm=i{E) ( 



(12) 



N{E) \2meJ as/E' 

where a = BD is the area of one period of the billiard. 

We see that Eq. (10) is precisely the first term in the Weyl series for the integrated density 
of states (see for example [12], Sect. 7). Equation (12) indicates that i) for fixed E and a the 
portion of m = 1 states, relative to the total number, is larger the narrower the billiard is and ii) the 
convergence of Nm=i{E)/N{E) to zero is rather slow. These estimates are in excellent agreement 
[22] with detailed numerical calculations for the case Nt = 1 and e <C 1 (e.g., e = 0.06). 

For the ease of a modulated boundary composed of many harmonics, the expression for Njn=i is 
still valid. However, N{E) is not simply given by (10) since now the perimeter 7 of the boundary, 
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as well as its curvature, contribute significantly to A''(£') as [30], 




The coefficients C2r+i, r > 1, in the sum of the Weyl series (13) depend on the curvature (and its 
derivative) of the modulated boundary (see Table 3, p. 545 of [30]). As is shown in [30], the higher 
the index r of the Wcyl coefficient C2r +i, the more complex the expansions for Cr arc, involving 
higher and higher powers of the curvature and its derivatives. Clearly, the larger the number of 
harmonics Nt is the more important the perimeter and curvature terms become. 

3 Localization Lengths of Eigenstates versus Profile's Com- 
plexity 

A detailed analysis of the eigenstates for the billiard studied here, with rippled or rough boundaries, 

shows that there arc three distinct types of eigenstates, namely, weakly, moderately, and strongly 
perturbed. To get a panorama of how these different types of states appear as a function of the 
energy, one may employ one of various quantities that give a measure of the effective number of 
unperturbed states that form the eigenstates of the perturbed system. One such measure is the 
"entropy localization length" [21] defined by Ir = exp [— (H — Hgqb)] ~ 2.08exp(— H). Here H is 
the Shannon entropy of an eigenstate in a given basis, K = ^^""^ | |^ In | C" and '^goe 
is the entropy of a completely chaotic state which is characterized by gaussian fluctuations (for 
Lmax oo) of all Components C" with the same variance <] 1^>= l/L^nax- Fig- 3 shows 
Ih as a function of the level number a [equivalently, as a function of the integral density of states 
N{E)] for the cases of N^i' = 2. 4, 6. and 100 harmonics, respectively. In all cases wc see that with 
the level number a, Ih increases on average. This is expected since with the increase of energy 
the eigenstates are expected to be more ergodic. Also, these figures demonstrate two other features 
which are only proper of quasi ID and 2D billiards with small amplitude of ripplcness or roughness 
on the boundaries: i) the increase on the number of states with minimal Ih as the complexity of 
the profile increases and ii) the strong and rapid fluctuations of Ih as a function of a. The first 
feature is easy to understand if we recall that the expression (11) for Nm=i is valid regardless of the 
complexity of the boundary while N{E) does depend on the complexity. In fact, as Nt increases 
the second term in (13) becomes more influential and hence the ratio Nm=i/N{E) increases with 
Nt. 

For any strongly localized state it is possible to determine the pair (n, m) of the unperturbed 
state that determines the main component of such localized state. Fig. 4 shows two sets of strongly 
localized states that have been identified with the unperturbed states (n, m) = (—18, 1) and {n = 
— 18,m = 2). The first row of states corresponds to the m = 1 unperturbed states and the second 
row to the m = 2 states. Note that the level number a corresponding to these states shifts down 
gradually as Nt increases, in agreement with the observation made above that the ratio Nrn=i/N{E) 
increases as Nt increases. A quick glance of figures 4(al-cl) also gives the impression that all m = 1 
states are qualitatively the same. However, some striking differences between them come to light 
when viewed in the configuration representation. Fig. 5 shows the same sequence of states as of Fig. 

4 but in the configuration representation. In this figure we see the one maximum and two maxima 
of the wave functions which are characteristic of the m = 1 and m = 2 states, however note that the 
wave function for the case Nt = 100 has become drastically asymmetric with respect to y = d/2; 
it is "repelled" by the rough boundary. We shall analyze this unexpected phenomenon in the next 
section. 

The second feature revealed by the plots of Ih as a fimction of a (Fig. 3) namely, the strong and 
rapid fiuctuations, reflects the coexistence of localized, sparse and extended states, characterized 
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Figure 3: Localization measure Ih for the first 1000 ordered in energy eigenfunctions for (a) Nt = 2, 
(b) Nt = 4, (c) Nt = 6, and (d) Nt = 100. 
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Figure 4: Localized eigenfunctions in the energy representation for (al) Nt = 2 (a = 987), (bl) 
NT = 6{a = 983), (el) Nt = 100 {a = 391), (a2) Nt = 2 {a = 990), (b2) Nt = 6 {a = 986), and 
(c2) Nt = 100 (a = 396). 
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FIGURE 5 

(available in gif format) 



Figure 5: Eigenfunctions of Fig. 4 in the configuration representation | 'ip°'{x,y;k) p; (af) Nt = 2 
(a = 987), (bl) iVT = 6 (a = 983), (cl) Nt = 100 (a = 391), (a2) TV^ = 2 (a = 990), (b2) Nt = 6 
{a = 986), and (c2) Nt = 100 (a = 396). 




Figure 6: Typical consecutive eigenfunctions in the energy representation for the case Nt = 100. 
(a) a = 390, (b) a = 391, and (c) a = 392. 



respectively, by small, medium, and large values of Ih- It was argued in the previous section 
that the number of unperturbed states conforming a given exact state depends on the value of the 
transversal mode m of the original unperturbed state. Since neighboring states in the energy-ordered 
unperturbed spectra generally have completely different values of the transversal number m, say, 
n = J^maxifn = 1 and n = l,m = M^ax^ it follows that the perturbed spectra should show rapid 
and strong fluctuations of Ir as a function of energy or a. This is expected to be independent of the 
complexity of the boundary. For the case Nt = 100, Fig. 6. shows a typical sequence of consecutive 
levels illustrating that if we arbitrarily pick up two neighboring states, knowledge of the structure 
of one does not allow us to predict the structure of the next. 

Additional information can be obtained by considering the reverse process, namely, the expansion 
of an unperturbed state in the basis of the perturbed system. In particular, we can study the "width" 
la or mean square root of an unperturbed state (j)^ defined asl^ = { Y^^^T I C'" |^ [a — nc(/)]^}^/^, 
where Uc = oi \ C" p gives the centroid of an unperturbed state in the perturbed basis. Fig. 7 
reports the plot of Icr as a function of I for Nt = 2, 4, 6, and 100. One can see a clear regularity in 
these plots. The lowest row corresponds to unperturbed states with m= 1 and m = 2. The next row 
correspond to imperturbed states with m = 2 and to = 3, and so on, see [21] for details. We see that 
as Nt increases the pattern of pairs persists for higher and higher values of the quantum number 
TO. This fact clearly reflects a regular character of the states with low values of to independently of 
how large the energy is. Note that the energy is not the unique semiclassical parameter, in fact we 
have two such parameters, to and n. 
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Figure 7: Localization measure for the first 1000 ordered in energy unperturbed states expanded 
in the basis of H for (a) Nt = 2, (b) Nt = 4, (c) Nt = 6, and (d) Nt = 100. 



4 Structure of Energy-Localized States for Disordered Pro- 
files 

Here we concentrate on the nature of eigenstates for the case of a rough boundary, Nt = 100. In Fig. 
8(al-b2) we present a typical pair of extremely energy-localized states together with their respective 
enlargements. Both are m = 1 states but with different values of n, namely {n,m) = (—18, 1), and 
(—28, 1). Close inspection of Figs. 8(a2-b2) shows that, unlike the case of profiles with small number 
of harmonics Nt, these states reveal other components of appreciable magnitude. A decoding of these 
components shows that the second most important component of each of these states corresponds 
to TO = 2, the next one to to = 3 and so on, all with the same value of n. That is, there is no 
mixing between different values of n. Moreover, as Figs. 8(a2-b2) suggest, their amplitude decay 
exponentially. We have found that, indeed, this is the case and hence one can write [see Eq. (3)], 



C?{k) = C^^ik) Sexp[-f3{m - 1)] S„ 



(14) 



where Ua corresponds to the energy-localized state. The parameters S = 0.8 and /? = 0.53 are 
obtained by fitting the data to the exponential dependence. It is important to remark that for a 
fixed Nt all energy-localized states give the same values of S and (3. Substituting Eq. (14) into Eq. 
(3) gives 



ipg^{u,v;k) = CSn-^/'^g-'^^^exp[i{k + nc)u] 

OO 

X ^ sin exp [-/3(m - 1)] 



where C arises to satisfy the orthonormality condition in curvilinear coordinates, 

p27V pi 



/ / dudvy^ip'j'Jip"^^ = 1. 
Jo Jo 



(15) 



(16) 
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Figure 8: Localized eigenfunctions in the basis representation (upper part) and an enlargement of 
the main components of them (lower part) for Nt = 100; (a) a = 398 and (b) a = 968. The 
main component corresponds to the unperturbed state characterized by (a) {n,m) = (18, 1) and (b) 
(n, to) = (—28, 1). The components shown in (a2-b2) correspond to unperturbed states with to =1, 
2, 3, 4, 5, 6, 7, and 8, respectively (from left to right). 



It follows that 



and 



^"V'/' ^ exp [-2/?(to - 1) 



^sin 



.m=l 



Km=l 



exp [— /3(to — 1)] 



(17) 



Cc(y)l' = TT-M £exp[-2/?(TO-l)] 



sin (J^) exp [-/3(to - 1)] 



_m— 1 



1 - exp(2/3) si-a^iny/d) 

A-K cos(7rj//d) — cosh(/3) 



(18) 



This latter expression gives the average shape of the localized eigenstates in configuration repre- 
sentation, projected onto the y— | V'" P plane. Also note that %l>f^^ does not depend on the parameter 
S, see Eq. (14). To obtain (18) we have made the approximation v ~ y and g ~ 1 since e = 0.06 <C 1. 
In Fig. 9 we plot Eq. (18) together with the numerical data for the state a = 391, see also Figs. 
4(cl) and 5(cl). The excellent agreement indicates that for extremely localized states with to, = 1 
there is no need to diagonalize the whole matrix since there is no mixing between different values of 
n. In order to obtain a good approximation of these eigenstates one only needs to diagonalize the 
block of size {M^ax x ^max) corresponding to the quantum number n^, which is a good quantum 
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Figure 9: Density | ip°'{x,y;k) p (points) for a = 391 and prediction (18) for strongly focalized 
eigenstates for Nt = 100: | tpgj{y) f (white line) [with /3 = 0.53]. Inset: | V'foc(y) 1^' ^q. (18), for 
(a) iVr = 1, /? = 0, (b) Nt = 25, P = 1, (c) Nt = 50, /3 = 1.3, (d) Nt = 75, P = 1.67, and (e) 
iVr ^ 00, /3 = 2.08. 



Figure 10: Localized eigenfunctions for Nt = 100 in the basis and configuration representation; (a) 
a = 976 and (b) a = 990. 

number in the presence of perturbation. In the inset of Fig. 9 we show global shapes of the m = 1 

states for different values of Nt- Wc have checked that these curves fit the m ~ \ states very well. 
Since the parameters S and /? are the same for all m = 1 states for a fixed Nt one can find S and 
/3 for some low energy state (so that the matrix to be diagonalized is small) and with this one can 
find the global shape of all m = \ states. Conversely, knowing the shape of a to = 1 state within 
any energy range, one may assess by looking at the amount of repulsion, the degree of complexity, 
presumably unknown, of the boundaries. 

Thus, the states most localized in energy, correspond to the smallest values of to of the originally 
unperturbed system. One can expect the same type of behavior, namely, the repulsion of the wave 
function away from the rough boundary for m = 2, see Fig. 10. However the main components of 
the wave function belong to different combinations of n and m. Nonetheless, their amplitudes decay 
roughly exponentially. One may treat the repulsion effect as the inability of the wave function to 
resolve the fine details of the rough profile and hence to satisfy the Dirichlet boundary conditions. As 
the value of m of the originally unperturbed state becomes larger, details of the profile can be resolved 
better and consequently the probability density should become appreciable in the neighborhood of 
the rough boundary. In the next section we shall see that this is indeed what happens. 
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FIGURE 11 

(available in gif format) 



Figure 11: Extended eigenfunctions in the basis representation (left) and in the configuration 
representation \ip"{x,y;k)\'^ (right); a = 986 for Nt = 1, a = 582 for Nt = 50, a = 407 for 
Nt = 100. 

5 Structure of the Energy- Extended States: Transition from 
Smooth to Rough Boundaries 

Here we study the nature of those eigenstates of the billiard with large values of localization lengths. 
These are selected from plots of entropy localization lengths, like those of Fig. 3, for a given number 
of harmonics. In Fig. 11 we present a typical set of this type of states in energy representation 
for Nt = 1, 50, and 100. All of these states were chosen to be of approximately the same energy, 
namely around the energy of the unperturbed level I = 980. Note that their level number a 
decrease as Nt increases because, as we showed in Section 2, the level number N{E) = a grows 
slower as a function of energy as the complexity of the profile increases. For smooth profiles with 
Nt ^ 1 the amplitude of its components decay roughly symmetrically around I = 980. However, 
as Nt increases the states spread more towards the high-lying unperturbed states with not a single 
predominant component. 

The strong asymmetry can be understood as follows. As the perturbation increases, more and 
more components arc needed to construct the exact states of the modulated billiard. Since the 
unperturbed levels are bounded from below, higher-lying unperturbed states are needed to form an 
exact state. Another argument rests on the classical-quantum correspondence of the structure of 
eigenstates [20, 21, 22, 31]. The structure of an eigenstatc a is given by the average W{E'^ \ E") over 
a number of states around the state a. In [21] it is shown how to construct the classical counterpart 
of the eigenstates (CCSE) for non-integrable billiards. Briefly, in analogy with the quantum case 
where we project the cigcnstate of the full Hamiltonian H on the basis of H'^, classically we project 
the dynamics of H onto Operationally, this is done by substituting the solutions of H into 
H^. Specifically, we write the Hamiltonian in the curvilinear coordinates {u,v). Thus the classical 
Hamiltonian for our billiard can be written as if = + eV, where, 

H' = ^{P^+P^) (19) 



and 

V = -- 



2me 



2m,_ 



(20) 



By substituting the solutions of H into H^, becomes a function of time, varying within an 
upper iJyp and a lower -fffo^^, limit. The energy range spanned by the projection of the dynamics of 
the full Hamiltonian onto if" is /S.E = H^^ — H^^^, referred to as the width of the "energy shell". 
Since the dynamics of H is fully chaotic, it is enough to substitute a single trajectory of H into 
and to let it run for a sufficiently long time to span the whole available phase space. From this, one 
obtains the distribution of energy values which is the CCSE in energy space [21]. 

In Fig. 12 we show the structure of some eigenstates and their classical counterpart for Nt = 2 
and 10. Quantum oscillations around the classical curve are found to be quite strong because the 
average generally includes extended and localized or sparse states [21]. Moreover, in average, the 
range of the energy shell and the shape of the CCSE correspond quite well to the global shape of 
eigenstates. 

Below we analyze the particular peak structure of these curves but before let us discuss the 
dynamical origin of the observed increase of H^^ as a function of Nt- We first note that the new 
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(a) 




Figure 12: Structure W{Ei \ E") of eigenstates in the energy representation (broken curve) and the 
CCSE (thick smooth curve), as a function of the normahzed energy /E"^. The average is taken 
over the interval (a) 2200 < a < 2400 for Nt = 2, and (b) 2200 < a < 2320 for Nt = 10. 



canonical momenta are given by 



and 



P^ = {l + eOPy, 



(21) 

(22) 



where Px and Py are both within the range [—y^2m^,^/2rn^E], with E being the energy of the 
total Hamiltonian. Prom Eqs. (21) and (22) we see that Pv is of the same order as Py but P„ can 
be much larger than P^ since the derivative = — J2n^ ^-^n sin(iVx) can be very large, roughly 
proportional to Nt for large Nt- 

The upper limit H^p of the energy shell then ought to grow roughly proportional to N^. Inspection 
of plots of Fig. 13 shows that indeed this is the case. For example, for Nt = 50 and Nt = 25, the 
ratio Hup,bo/ Hup,2b = (50/25)^ = 4, which agrees well with the limits found numerically. 

Now let us study the form of eigenstates as a function of Nt- For the case Nt = 2 in Fig. 12(a) 
we see that the CCSE has four peaks in a rough agreement with the global shape of eigenstates. 
For Nt = 1 in [21] it was shown that both classical and quantum structures have three peaks, a 
central one and one on each side of the central peak. It was found that the main contribution to the 
central peak comes from grazing trajectories (bouncing balls), while the right (left) peak is due to 
trajectories dwelling in the neighborhood of the period one stable (unstable) fixed point. The profile 
^(u) = Ai cos(u) + A2 cos(2u) for the case Nt = 2 gives rise to three periodic orbits of period one; 
two stable ones at {Pu,u) = (0,0) and (0,7r) and one imstablc at {Pu,u) = (0, cos^-'^(— Ai/4y42)). 
At these points, the values of the profile function arc: ^(0) = Ai + A2,^{n) = —Ai + A2, and 
^(cos~i(-Ai/4^2)) = -{Al/9,A2) - A2. Substituting P„ = and = into (19-20) one obtains 
P^ = 2m,,,H{l + e^) and iJ" = (1 + e!;)^H, where ^ stands for evaluated at the fixed period 
one points. Inserting the randomly generated values of Ai and A2 used in the calculations, we 
obtain i?" — 0.899, 1.0681, and 1.1409, corresponding to the periodic points at u = 0, tt, and 
cos~^(— A1/4A2), respectively. The two peaks on the right of the highest peak of the CCSE of Fig. 
12(a) correspond to the stable period one orbits and the left peak to the unstable period one orbit. 
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Figure 13: CCSE for Nt = 10, 25, and 50. 



Clearly, the number and location in energy of the peaks depend on the coefficients An- Then by 
looking at the structure of the eigenstates, one can infer about the complexity of the profile. As the 
number of harmonics determining the profile increases, AiJ grows proportional to and all the 
peaks produced by periodic orbits of period one, which are located very close to the central peak at 
= H, seem to disappear in favor of a wider peak around = H. 

Thus, the CCSE allows us to understand characteristic features of the eigenstates in terms of the 
underlying classical dynamics. 



6 Energy Band-Structure: Localized versus Extended States 



We address here the question of how the energy-localized and energy-extended states show up in 
the energy band structure of the billiard as the complexity of the profile increases from smooth to 
rough. 

We have seen that the energy localized states v; k) are well described by Eq. (15). Since 

they are very close to unperturbed states, their normalized energy E"' / E* is given approximately 

by 

if \ 



E"/E* 

and the group velocity (the slope of the bands) is 



IT 



d{E°'/E* 
dk 



-^{k + no) = 8(fc + no 



(23) 



(24) 



where E'^ ^ El^^{k) and E* = {J)f /2me){-K/df. 

In the last equation we substituted d = 27t, the value used throughout this paper. The group 
velocity Vg as given by Eq. (24) should be very close to the exact one regardless of the value of Nt- 
To check our expectation let us analyze Fig. 14 which shows a segment of the energy-band structure 
for the case Nt = 100 in the range 0.05 < k < 0.15, and 1280 < E"/E* < 1370. Note that there 
are mostly only two types of bands, flat bands and straight sloped bands, everywhere except in the 
neighborhood of avoided crossings. In this figure we indicate the position of some representative 
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eigenstates at A; = 0.1 which are shown in their energy and configuration representation in Fig. 
15. Inspection of these two figures reveal that: i) eigenstates lying on horizontal (inclined) bands 
correspond to energy-extended (localized) states, and ii) the slopes computed directly from the 
graphs agree very well with the group velocity given by Eq. (24). Therefore, the inclined bands are 
actually parabolic bands but because the range in k is small, they appear in this figure as straight 
sloped lines. The agreement is quantitative, for example, the slope of the bands at fc = 0.1 for states 
391 and 408 is ±245.5. The use of Eq. (24) with k = 0.1 gives = 18.09, with the integer part 18 
being exactly what is found by identifying the (n, m) pair of the unperturbed state that coincides 
with the main component of the energy-localized state (a = 391 or 408). Such a good agreement 
is surprising for the case of the state a = 391 since it lies near an avoided crossing, where one may 
anticipate strong mixing of many unperturbed components. We have examined the structure of the 
states in the neighborhood of the avoided crossing as we vary k (see inset of Fig. 14) and observed 
a gradual exchange of identities between states. Specifically, for k = 0.1022 the state a = 390 (391) 
is extended (localized) in energy, and as k increases the state a = 391 acquires more and more 
components till at A; = 0.1034 it becomes extended throughout the whole energy shell, looking very 
much like the state a = 390 at = 0.102. Thus, the states retain their character (extended or 
localized) throughout the Brillouin zone except in a small neighborhood of the avoided crossings. 




0.05 0.1 0.15 

k 

Figure 14: Energy band structure for Nt = 100, and (a/27r, d/27r) = (0.06,1.0). Inset: enlarged 
avoided crossing marked as a black rectangle at the right of the states a = 390 and 391. 

It should be stressed that for any high energy one can find eigenstates with a very different 
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FIGURE 15 

(available in gif format) 



Figure 15: The eigenstates a = 390, 391, 407, and 408, marked in Fig. 14, in the basis representation 
(left) and in the configuration representation (right). 
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Figure 16: Entropy localization length Ih and group velocity \Vg\ for the eigenstates in the range of 
energy of Fig. 14, a = [385,414] {Nt = 100, k = 0.1). 



structure. Note that the energy-extended states of Fig. 15 are strongly localized in the coordinate 
X within one period of the billiard. We found that in the case of rough profiles, about a third of 
energy-extended states become localized in configuration representation. This localization effect is 
in agreement with the structure of the corresponding energy band spectra: since these states have 

~ their wave functions are not allowed to propagate through the billiard. Also from Fig. 11 
one can clearly see that the transition from smooth to rough boundary in the billiard leads to a 
transition from extended to localized in configuration representation of the energy-extended states. 

Since states located on flat (parabolic) bands correspond to energy-extended (-localized) states, 
their localization lengths are expected to be roughly proportional to | |~^. A detailed comparison 
shows that while this inverse proportionality is qualitatively correct, the group velocity docs not 
distinguish between fully and partially energy-extended states, they both have | Vg \ close to zero. 
This is illustrated in Fig. 16. Nevertheless, it is clear that the group velocity gives useful information 
about the nature of a given state and, besides, it is the band structure which is usually known, and 
not the wavefunction itself. Thus, a graph of as a function of a, for an arbitrary value of k within 
the Brillouin zone can be very helpful in understanding the global structure of the energy band and 
the type of corresponding eigenstates, localized or extended in energy space. 

Fig. 17 shows Ih and | | as a function of a for Nt — 1, 25, 50, 75, and 100. We consider first the 
case Nt = 100, noting the dominance of states with zero or near zero group velocities (equivalently, 
maximal values of Ih)- This is consistent with Fig. 14 but Fig. 17 also strongly suggests that both 
types of bands (fiat and parabolic) and, correspondingly, the both types of eigenstates (extended 
and localized in energy) appear at all energy ranges, with the exception of the lowest levels where the 
perturbative regime holds. Further, the sequence of graphs in Fig. 17 reveal that the distribution 
of slopes or entropy localization lengths changes gradually as the complexity of the profile changes. 
Ona can compare the plots of |V^| as a function of a of Fig. 17 for Nt = 1 and 100. The former 
shows the no-proliferation of flat bands, in fact, the values of | | appear evenly distributed: there 
are more intermediate than maximal or minimal values of | Vg | . It would imply fluctuations of the 
bands yielding an spaghetti-like structure. This is indeed the case which is qualitatively different 
from that of the rough profile discussed above. The rest plots of Fig. 17 imply that there is a 
gradual transition in the energy band structure as Nt increases and this in turn means that as 
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Figure 17: Entropy localization length Ih (left) and group velocity \Vg\ (right) for the first 1000 
eigenstates for Nt = 1, 25, 50, 75, and 100; k = 0.1. 



the underlying classical dynamics goes from chaotic to random the eigenstates with intermediate 
localization lengths tend to disappear in favor of either extremely extended (in energy) or extremely 
localized (in energy) states. Hence, knowledge of the group velocity gives information not only about 
the type of energy band spectra but also about the underlying classical dynamics, in its transition 
from chaos to disorder. 



7 Conclusions 

We have studied the transition from chaos to disorder in periodic billiards, focusing on the properties 

of eigenstates and energy band spectra. We did this by considering a prototype 2D billiard, defined 
by two parallel walls, one wall is flat and the other is defined by a sum of Nt harmonics with 
randomly distributed amplitudes. The larger the sum the more complex the profile, and for a large 
Nt the profile becomes rough (disordered). We remark that we considered the ease where the 
root mean square of the amplitude a is small compared to the the average width d of the billiard, 
{a/d-^l). 

We found that independent of the number of harmonics, and for any energy range, there are 
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extremely localized states co-existing with sparse and extended states (in the unperturbed energy 
basis). These "energy- localized states" originate from unperturbed states with the smallest values 
of the transversal mode m = 1, 2, 3. Being very similar to traveling waves, these states have their 
maximum support along one (m = 1) or two (m = 2) horizontal strips. Hence they are identified as 
bouncing ball states. As the complexity of the profile increases, these strips are gradually repelled 
towards the flat wall. Hence, knowledge of the amount of repulsion gives an indication of how smooth 
or complex the modulated profile is. 

We showed that the number of the energy-localized states relative to the total number of states 
up to some energy E is proportional to E~^/'^, hence the approach to ergodicity is rather slow. 
Further, for small Nt, these energy-localized states coexist in roughly equal numbers with sparse 
and extended states but as the complexity of the profile increases the number of sparse states 
decreases in favor of ergodic states. We found that in the case of rough profiles, about a third of 
energy-extended states become localized within one period of the billiard in configuration space. 

While analyzing the energy band structure, we showed that the "dynamics" of the levels for the 
case of rough boundaries is quite simple: they follow either horizontal or parabolic trajectories in 
/c-space, except very near the avoided crossings. States lying on horizontal (parabolic) bands are 
energy-extended (localized) states. In contrast, for smooth boundaries Nt ~ 1, the bands form the 
typical spaghetti-like structure of chaotic systems. Thus, the energy band structure reflects clearly 
the degree of complexity of the boundaries and, in turn, the transition from chaos to disorder. It was 
also shown that the group velocity Vg = dE/dk is a very useful quantity for understanding dynainical 
features of the system. Specifically, plots of as a function of a, where a labels the energy-eigenstate 
of the perturbed Hamiltonian, show distinct patterns for smooth and rough boundaries. The utility 
lies in that it is only necessary to calculate as a function of a for any value of fc, to predict 
the type of energy band structure and also the type of underlying classical dynamics (chaotic or 
disordered). 
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